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We investigate effects of optical lattice potential in one- and two-dimensional two-component 
trapped Fermi gases with population imbalances. Using the exact diagonalization and the density 
matrix renormalization group methods complementarily, we calculate the atom density profile from 
the ground state many-body wavefunction as a function of attractive interaction strength for various 
population imbalances. The numerical results reveal that although a phase separation between the 
superfluid core and the shell cloud of excess atoms occurs as observed in experiments without the 
optical lattice, the population imbalance generally remains in the core region in contrast to the non- 
lattice cases. The essence of the numerical results in a strong attractive regime can be explained by 
an effective model composed of Cooper pairs and excess major fermions. 

PACS numbers: 03.75.Ss, 71.10.Fd, 74.81.-g, 74.25.Jb 



I. INTRODUCTION 

Recently, two-component fermion systems with popu- 
lation imbalance have attracted much attention in vari- 
ous research fields, such as cold atoms, superconductors, 
and QCD Q. In the 1960's, effects of the population 
imbalance have been theoretically investigated in the su- 
perconductivity literature [3, d, [1, H| • Sarma considered 
the stability of the gap phase (Sarma state) [|J and Liu 
and Wilczek revisited it with a new picture (interior gap 
phase) PL Fulde and Ferrell [H, and Larkin and Ov- 
chinikov [J] predicted the so-called FFLO state, where 
the superconducting order parameter is spatially modu- 
lated. Very recently, some evidence of the FFLO state 
has been reported in a heavy fermion superconductor 
CeCoIn 5 0. 

The population imbalance has been also extensively 
studied in ultra-cold Fermi gases @, S, B 0, HH, O ■ The 
advantage of using atom gases is that one can widely tune 
some physical parameters, such as the interaction and the 
population imbalance. In a trapped two-component 6 Li 
Fermi gas @, H, [1, EU E3 > a phase separation between a 
superfluid core region and a surrounding unpaired gas of 
excess atoms was observed. Recently, the absence of the 
population imbalance inside the superfluid core region 
was reported in 0, [ll[ . 

So far, the superfluid Fermi gas with population imbal- 
ance has been experimentally studied without an optical 
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lattice 0, H, OH [IH H2I. In th is p aper, we investi- 
gate effects of the lattice [13l Ibl . ITa . fig on the pop- 
ulation imbalanced systems. Very recently, there have 
been several reports on this issue, one of which theoreti- 
cally claims that FFLO is more stabilized than the non- 
optical lattice case , and others of which numerically 
show that FFLO is observable in the one-dimensional 
attractive Hubbard model by using the density ma- 
trix renormalization Group (DMRG) method [HI, I20I ]. 
Our approach in this paper differs from these works 
in using the exact diag onalization method mainly and 
the DMRG method [l9l. |2C| complementarily in order to 
check whether the exact diagonalization results are small 
size effects or not. In addition, we focus on the density 
profile of majority (nj) and minority components (rij) 
at T = in a one- and two-dimensional lattice trapped 
systems instead of the stability of the exotic superfluid- 
ity FFLO. The reason is that the atom density profile 
is the most convenient and clear observable for experi- 
ments. Moreover, the FFLO requires at least the popu- 
lation imbalance in the atom density profile. Namely, it 
is crucial point for FFLO whether the population imbal- 
ance survives even in the presence of pair instability at 
T = ground state or not. Thus, we widely measure the 
atom density profiles in various system sizes, imbalance, 
interaction, dimensionality (1D-2D), and traps, prior to 
an examination of the FFLO stability. Consequently, we 
point out in a strongly-attractive interaction regime that 
correlation effects due to the presence of lattice result in 
a core phase with a typical population imbalance, which 
is very different from non-lattice cases [1, [ll[ . 

The contents of this paper are as follows. In Section 
HTl we summarize the experimental results of the atom 
density profiles in the imbalanced Fermi gases without 
the optical lattice and show our prediction on those of 
the imbalanced gases in the presence of the optical lat- 
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FIG. 1: Schematic profiles of the total atom density and the 
density difference between particles with pseudo-spin up and 
down, i.e., n-f — n±, which are measured along the long axis of 
the trap for polarized fermionic gases. The upper panels (a) 
and (b) display the experimental (non-lattice system) results, 
while the lower panels (c) and (d) predict the profiles in the 
presence of the optical lattice. The parabolic lines schemat- 
ically represent the trap potential shape. In the absence of 
the optical lattice, the experiments revealed that the phase 
separation between the paired (superfluid) gas and the com- 
pletely polarized (excess major atoms only) periphery occurs 
as shown in the upper panels (a) and (b) . On the other hand, 
the lower panels (c) and (d) show a typical result of numerical 
calculations in the presence of the optical lattice. A big differ- 
ence between the upper and the lower panels is that the excess 
species can stay even on the superfluid core in the presence 
of the optical lattice. 

tice. Section IIIII is devoted to explain the model and 
the numerical method. In Section IIV1 numerical results 
by the exact diagonalization and DMRG are presented, 
theoretical analysis using an effective model derived in 
strong coupling limit is given, and strong trap effects are 
displayed as an exceptional case. Section [V] shows the 
results in two dimensional (2-D) systems, and Section [VTl 
discusses temperature effects on all the present results. 
Finally, summary and conclusion are given. 



II. PARTICLE DENSITY PROFILES IN 
POLARIZED ATOMIC FERMI GASES 

In this section, we summarize the experimental results 
[U of the trapped imbalanced Fermi gases without 
the optical lattice and predict how their results change 
by adding the optical lattice. Figure Q] (a) and (b) are 
schematic pictures of atomic density profiles obtained in 
experiments without the optical lattice, and (c) and (d) 
are our prediction of those in the presence of the optical 
lattice. At first, let us focus on Fig.[T](a) and (b). The ex- 
periments revealed that the polarized gases show a phase 
separation between the core and the periphery below the 
superfluid transition temperature. The phase separation 
picture was clearly confirmed by monitoring the atom 
density profile for each of two Fermi- atom species. In- 
side the core, the profile shows an atom density ratio, 
n.f(x) : ni(x) = 1 : 1 as shown in Fig.[T](b) [§, llH, while 
a completely-polarized phase appears in the periphery, 
where the atom density ratio, n^(x) : n[{x) = 1 : in 



a case that the majority species is Fermi atom with the 
pseudo-spin \. 

On the other hand, our calculations in the presence of 
the optical lattice reveal that although the phase separa- 
tion occurs like the experimental results without the opti- 
cal lattice, the atom density profile shows a big difference 
compared to the non-lattice cases. We find in the pres- 
ence of lattice that the population imbalance generally 
remains even in the superfluid core region, i.e., the atom 
density ratio inside the core, n^{i) : ^ 1 : 1. More- 
over, we find in sufficiently strong pairing-interaction 
that ni(i) : n±(i) ~ 2 : 1 [21j inside the core as shown 
in Fig. Q] (d). This result is compared with an effec- 
tive model theoretically derived in the strong limit of 
the attractive interaction. In this paper, we claim that 
correlation effects on the lattice bring about the imbal- 
anced core phase. The main correlation effects on the 
lattice originate from a repulsive interaction between two 
Cooper pairs [12] and a site-exchange energy gain be- 
tween a Cooper pair and an excess fermion as shown in 
IIV1 The effective model successfully explains why the 
imbalanced core emerges. 



III. MODEL AND NUMERICAL 
CALCULATION METHOD 

Let us show the model Hamiltonian and the numer- 
ical calculation method to obtain the ground state of 
the model. In this paper, we consider two-component 
trapped Fermi gases inside optical lattices with N lattice 
sites. When the optical lattice potential is sufficiently 
strong, the system can be well described by the Hubbard 
model, the one-dimensional (1-D) Hamiltonian of which 
is given as (HHl], 

^Hubbard = -t ^ i a ]<j a i<? + H ' C + U ^2 Ui 1 ni l 




where a ia is the creation operator of a Fermi particle 
with pseudo-spin a —\ or J, at i-th site, and U (< 0) 
is a pairing interaction. The summation is taken 

over the nearest-neighbor sites, and the hopping matrix 
element t describes the hopping energy of atoms between 
the nearest-neighbor sites. The last term in Eq. ([1]) gives 
a harmonic trap potential, in which V corresponds to the 
potential height at the edge of the lattice. At the edge, 
the open-boundary condition is imposed. If the atomic 
density profile drops down at the edge, then the bound- 
ary condition does not almost affect the result. The pa- 
rameter V and the total number of Fermi atoms Np are 
selected to sufficiently reduce the atom density at the 
edges. 

In this paper, we examine various population imbal- 
ances with keeping Nf > N± (Np = ATf + N±), where JVj 
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and Ni are the total number of Fermi atoms with pseudo- 
spin | and |, respectively. In 1-D cases, we numerically 
diagonalize the Hamiltonian, Eq. {1} to calculate the den- 
sity profile at T = 0. Although this approach gives us ex- 
act results, the accessible system size is severely limited. 
To compensate this disadvantage and confirm whether 
the exact diagonalization results are small size effects or 
not, we employ the DMRG method. The DMRG guar- 
antees a high accurate result as long as the trap potential 
is not so steep. We confirm that a difference between the 
eigen-values obtained by both methods is below 1CP 6 in 
the lattice size N = 20 with V/t = 1. Since the DMRG is 
mainly used for much less steep traps (V/t is fixed to be 1 
except for Section VI D, which examines a trap strength 
dependence), the DMRG results sufficiently keep high 
accuracy. 

On the other hand, in two-dimensional (2-D) square 
lattice systems, we use the exact diagonalization solely, 
because the method is now the most accurate and reli- 
able for 2-D finite systems. Here, we just mention that 
the direct extension of the DMRG to 2-D ladder cases, 
which keeps high accuracy like 1-D cases is now under 
development [25|, [2fj| . For the exact diagonalization and 
the DMRG, we use SX-6 (1-node 4 CPU's system) and 
Altix 3700Bx2 (scalar parallel machine) in JAEA and the 
Earth Simulator [27|]. In the exact diagonalization, the 
problems with the lattice size N — 16 are calculated by 
SX-6, while those with above N = 20 the Earth Sim- 
ulator. In the latter case, we need parallelization and 
high-performance computing techniques. See Ref. [27j 
for computational technical issues on massively parallel 
supercomputers. On the other hand, all DMRG calcula- 
tions are made on Altix 3700Bx2. 



IV. NUMERICAL RESULTS FOR 1-D 
IMBALANCED GASES 

A. Exact Diagonalization Studies 

Let us present exact diagonalization results for the im- 
balanced 1-D attractive Hubbard model with the trap 
potential. Figure [2] shows U / 1 dependences of the cal- 
culated atom-density profile of the ground state, n(i) = 
n-\(i) + n[(i) for three population imbalances, whose ra- 
tios of major to minor species are (a) 5:3, (b) 6:2, and 
(c) 7:1, respectively. For these three cases, the lattice 
size N = 16, the total number of fermions TVp = 8, and 
V/t = 1 . The filling is just the quarter one if the trap po- 
tential is absent. Although the depression at the edges is 
not completely zero in the presence of the trap (V/t = 1), 
we confirm that there are no edge effects because of no 
significant change of profiles on the lattice size extension 
(to N — 18 and 20) with the same fermion number (see 
Fig. [2](d)). 

In each panel, we find a step-like structure (see arrow 
in each figure) in each dome like profile. This structure 
becomes more pronounced as one increases the amplitude 




FIG. 2: The exact diagonalization results for U /t dependences 
( — 15 < U/t < 0) of particle density profiles, n(i) for three im- 
balance cases, whose ratios of major to minor fermion species 
are (a) 5:3, (b) 6:2, and (c) 7:1, respectively. In these three 
cases, Np = 8, N = 16 and V/t = 1. The panel (d) is a 
highlight of n(i) at U/t — —10 for the imbalance from 4:4 to 
7:1 in TV = 20. Nf and V/t are the same as (a-c). 



strength of the attractive interaction |J7|. The structure 
is also observable in Fig. 2(d) with a larger size (N = 20), 
in which Np and V/t is the same as (a-c), and U/t is 
fixed to be —10. 

The step-like structures seen in Fig. [2] are found to 
be associated with the phase separation between a core 
phase whose main component is pair and a shell one com- 
posed of only un-paired excess atoms from other observ- 
ables as shown in Fig. [3] Figure [3] (a), (b), (c) and (d), 
respectively, show U/t dependent profiles of the single- 
occupation density np(i), the double-occupation density 
ne(«), the density subtraction Jij(i) — ni(i), and the on- 
site pair amplitude A(i) 0, HU defined by 

= ^(A T +1,N L + l|aIXjiV T , N t ) (2) 

where, | JV-f, JVj_) is the ground state wavefunction with 
(iVf,.iVjJ atoms. This function is a probe for the local 
condensate amplitude [28[ for finite systems. Here, we 
note that the probe gives non-zero values even at the 
interaction free (U/t — 0) as seen in [3] (d) and shows 
spatial modulations without sign changes especially in a 
relatively-weak interaction regime. The reason for the 
former one is that the function counts all the possible 
pair formation probabilities on a local site, i.e., it in- 
cludes even the pair formation due to the trap potential, 
and that for the latter one may be that the strong con- 
finement effect masks the sign change due to the FFLO. 
Although this latter interpretation requires further ex- 
aminations, we confirm at least the growth of the local 
condensate amplitude. See Ref. [l8| for the spatial pair 
correlation in larger core phases. n-p(i) and n&(i) are 
calculated by picking up the coefficient of the eigen-state 
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whose i-th site occupation is 1 and 2, respectively, from 
the ground state wave- function and summing up the am- 
plitude of the coefficients. Their explicit operator defi- 
nitions are given by np(i) =< n^(i) + Uj,(i) > — 
where ns(i) =< n^(i)ni(i) >. For the methodological 
details to calculate nF(i) and ub(j), see Ref. [I2|. Ac- 
cording to Ref. [22^ . as |J7/t| increases in the balanced 
case, nv{i) decreases, while 713(1) instead increases. This 
simple relationship means that the Cooper pair becomes 
more tightly-bounded one on a site with increasing |E//t|. 



From Fig. [3j we clearly find that 71b(«) and A(i) well 
develop only inside the "phase boundary" characterized 
by the step-like structure as seen in Fig. [2] while they 
completely diminish outside the boundary (shell region) 
(2^ |. This tendency becomes more pronounced with in- 
creasing the attractive interaction and somewhat satu- 
rated in the strong limit. This suggests that the pair 
density and the superfluid order (phase rigidity) locally 
develop inside the core. This does not contradict a con- 
sensus that the superfluid correlation grows as the most 
dominant one in f-D attractive Hubbard model [30j al- 
though any correlations algebraically decay in 1-D sys- 
tems. Thus, we call the core region "superfluid core" in 
the following. However, we note that there is a subtle 
issue whether the whole core really exhibits superfluidity 
in non-zero temperature due to its dimensionality and 
finiteness or not. The discussion originates from a prob- 
lem whether the phase rigidity grows and survives from 
an edge of core to another edge against fluctuations. This 
is beyond information obtained by the present analysis of 
the exact diagonalization and DMRG. Thus, we do not 
discuss the stability of the superfluidity but concentrate 
only on atom density profiles in the ground state in this 
paper. We believe that the future experiments and more 
advanced theoretical and numerical methods resolve the 
problem (e.g., see Ref. [l8[ for the static non-local corre- 
lation of the pair function inside the core phase). 



Let us turn back to atom density profile in the pres- 
ence of lattice, again. Fig. [3] (a) reveals that un-paired 
component stays even in the core region, where the 
ne(i) and A(i) grow. As a result, the density differ- 
ence 5n(i) = nf(i) — n-i(i) does not vanish in the core 
region as seen in Fig. [3] (c). Generally, we find that 
n-f(i) : ni(i) ~ 2 : f in the core region except for a very 
weak coupling regime (compare Fig.[^b) with Fig.[3]Jc)). 
This feature is clearly seen in Fig.4(a-c), which display 
slices of n-[{i") and nj,(i) profiles for U/t = —15, —5, and 
— 1, respectively. Inside this range, it is found that the ra- 
tio is approximately 2:1. Also, the ratio is seen for other 
imbalance cases (see e.g., Fig. 111(d)), although there is a 
small deviation from 2:1. This is quite different from the 
experimental result in the absence of the lattice, where 
one obtains Jif(i) : n±(i) ~ 1 : 1, or 5n(i) ~ in the core 
region @,|ll|. 




FIG. 3: The exact diagonalization results for U /t dependences 
( — 15 < U/t < 0) of profiles of (a) the single occupation den- 
sity riF(i), (b) the double occupation density riB(i), (c) the 
subtraction from major to minor species density n,j (i) — nj. (i), 
and (d) the on-site pair amplitude A(i). In all the cases, 
Nf — 8 (6 j, 2 I), and other conditions are the same as those 
of Fig.Q](a-c). 




FIG. 4: The exact diagonalization results of n-f(i) and ni(i) 
profiles for (a) U/t = —15, (b) —5, and (c) — 1. The numerical 
condition is the same as the case of Fig. 2(b), in which the 
imbalance ratio is 3 : 1. (d) rif(i) and nj_ (i) profiles for U/t = 
— 15 in the case of the imbalance ratio 7 : 1. 



B. DMRG Studies 

In order to confirm the generality of the phase sepa- 
ration and the core phase characterized by the approx- 
imate ratio "2 : 1," we adopt the DMRG method and 
study much larger systems (N = 80 and 160). Figured] 
(a) and (b) show U/t dependences (—15 < U/t < 0) of 
density profiles for major and minor species, respectively, 
and Fig. [5] (c) is their subtraction. In their results, the 
lattice size N — 80, Np — 32, and the prepaired ratio 
of different atom species is 3 : 1. It is found from Fig. [5] 
(a-b) that the atomic density ratio in the core region is 
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FIG. 5: The DMRG results of U/t dependences (-15 < U/t < 
0) for density profiles for (a) major species n-j- (i) , (b) minor 
ones n|(i), and (c) the subtraction from major to minor one 
n T (i) -n L (i). In these panels, N = 80, and N F = 32 (24 |,8 | 
). The panel (d) and (e) are two slices for U/t = —15 and —5 
for ri](i) and ny(i), and the panel (f) is a half left-side profile 
of the local density of n j ( i) and nj, (i) for TV = 160, Nf ~ 62 
(T 48, | 16), and U/t = -10. In all cases, V/t = 1. 



about "2 : 1" except for the weak coupling region. The 
Fig. [5] (d) and (e) are two slices (U/t = —15 and —5) of 
profiles of ri| and n\, which clearly exhibit the approxi- 
mate ratio "2 : 1". In addition, Fig. 0(f), which shows 
DMRG result (a half profile) of n T and n l for N = 160 
and Np = 64, is also a clear demonstration of the approx- 
imate ratio "2 : 1." In these cases, it is noted that the 
trap potential and the resulting confinement is relatively 
very weak compared to the exact diagonalization results. 
From the common results by both numerical methods, 
the approximate ratio rij : ~ 2 : 1 is widely expected 
to be observed in the strong coupling regime in the pres- 
ence of the optical lattice with harmonic traps. 



C. Effective Model and Discussion 

Both the exact diagonalization and the DMRG reveal 
that the excess un-paired atoms remain even inside the 
core region in contrast to the experimental results with- 
out the optical lattice. Especially, the ratio of majority 
to minority species is approximately "2:1" inside the core 
in the strong coupling regime. This result is considered 
to clearly originate from strong correlation effects due to 
the presence of lattice. Here, we derive an effective model 
in order to explain the findings. In the strong-coupling 
limit \U\ — > oo, one finds that the present lattice Fermi 
gas is effectively described by a gas mixture of unpaired 
free Fermi-atoms and dimer atoms [22|, [3l[ . The effective 



Hamiltonian has the form, (dropping the potential term) 

H Hubbard — * 2t 2 /U (n i+ i rij — Bj +1 Bi) 

i 

-^("i+Lt^l +H.C.) 

i 

'E^lu^+^Kt+H.C.). (3) 

i 

Here, B\ is a creation operator of a bound pair at the i- 
th site. In this effective Hamiltonian, the first term gives 
the nearest-neighbor repulsive interaction U c s = 2t 2 /U 
between Cooper-pair bosons and the energy gain of the 
molecular hopping — t e g = —2t 2 /U between the nearest- 
neighbor sites. The second term in Eq. stands for the 
kinetic energy of excess "["-spin atoms, and the third term 
describes an exchange energy between a Cooper pair and 
an excess atom when they sit on the nearest neighbor 
position each other. The third term is also derived by 
considering the hopping of the minor atom from an on- 
site Cooper pair to the nearest neighbor site, in which an 
excess major atom solely stays. This effective model is 
valid in |J7| — * oo, since there exists conversion process 
between the bosons and the excess fermions for finite neg- 
ative U/t. However, it is numerically found that the limit 
description (\U\ — > oo) becomes sufficiently good more 
than \U/t\ ~ 10, since the un-paired fermions almost dis- 
appear in the range. See Ref. [22j for U/t dependence of 
the total number of un-paired fermions (= Yli n F(i)) i n 
the balanced cases. 

In the effective model considering only the first term of 
Eq. ([3]), i.e., in the balanced cases, it is well-known that 
the pair density correlation function shows a large mag- 
nitude between every other sites. The reason is that the 
Cooper pair tends to occupy every other site to avoid 
energy loss by the nearest-neighbor repulsion between 
Cooper-pair bosons (see the first term of Eq. ©). As 
a result, such a correlation effect causes a static spatial 
modulation of the pair density in the presence of the trap, 
which breaks the translational symmetry. In other words, 
a degeneracy between two dominant states having the 
site density correlation of pairs as " — n-j-j, — — n-^ — 0—" 
and "—0 — n-|-| — — n^— " , where n-j^ is the site density 
of the Cooper pair (n^ < 1), is broken by the trap po- 
tential, and the static density modulation as represented 
by "— nil — n'li — nil — n"i" appears. This behavior in 
the balanced case was also confirmed by different authors 

[HIM Hi. 

On the other hand, in the imbalanced systems, both 
the first and third terms of Eq. ([3]) explain the present 
numerical results of the density profiles. The first term 
prefers two degenerate density correlations described 
above, and excess atoms can then find their positions 
in the nearest neighbor sites of the on-site pair. Fur- 
thermore, the third term stabilizes the two density- 
correlation profiles, in which excess majority atoms slip 
in like "— n^— n^— n^— n-\— " and "— n-\— n^— n-f— n^— ", 
since the energy 2t is gained by an anti-bonding exchange 
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60 

FIG. 6: The DMRG results of V/t dependence of the particle 
density profile for N = 60, N F = 40 (30 T, 10 |), and U/t = 
-10. 



between an on-site pair and its neighbor excess atom. As 
a result, one finds that the density difference in the core 
region is simply given by n-\ : ~ 2 : 1 in the large 
\U\ range. We emphasize that the presence of the lattice 
is responsible for the approximate ratio 2 : 1 while it is 
absent in a Fermi gas with no lattice potential. 



D. Strong Trap Effect 

In this subsection, we study a confinement force de- 
pendence on the atom density profile. Figure [6] shows 
V/t dependence on the density profile for the imbalance 
ratio TVj : N± — 3 : 1. The focus of interest is a case in 
which the strength of the trap potential is so strong that 
the confinement force overwhelms the correlation effects 
due to the presence of lattice, i.e., the repulsive interac- 
tion between Cooper pairs and the exchange energy gain 
between a Cooper pair and an excess atom as explained 
above. Such an extreme situation is observable above 
V/t ~ 35 as seen in Fig. [51 One finds above V/t ~ 35 
that un-paired atoms are excluded from a central part of 
the core region. Then, the density difference from major 
to minor species becomes Sn(i) ~ 0, namely, the ratio 
»if : 71^ ~ 1 : 1. This feature is equivalent with the re- 
cent experiments without lattices [a, El- However, such 
a situation may be non-realistic, since it requires very 
large V compared to the non-lattice one. Moreover, the 
numerical simulation result (Fig. [6| reveals that only a 
central point of the core region shows nj : rij = 1 : 1. 
This is also not experimentally observable in non-lattice 
cases. It is found even in strongly confined cases that the 
lattice presence is non-negligible. 



V. NUMERICAL RESULTS FOR 2-D 
IMBALANCED GASES 

In this section, we study 2-D imbalanced gases in the 
presence of the optical lattice. The lattice is a square one, 
in which only the nearest neighbor hopping is allowed. 
In this paper, we employ only the exact diagonalization. 
Thus, the total number of lattice sites N is limited to 25, 




(c) Bt(i)-»4.(f) 



FIG. 7: The exact diagonalization results of the particle den- 
sity profiles of (a) n-j- (i) , (b) and (c) nj(i) — n±(i) for 2-D 
attractive imbalanced Hubbard model with a 2-D harmonic 
trap potential. In this case, N = 25, Nf = 12 (9 f, 3 j), 
U/t = -15, and V/t = 1. 



and Np = 12, where the imbalance ratio n-f : n± = 3 : 1 
(9T,3|). 

Figure[7]shows a typical result of the calculated density 
profiles for n-\(i), ni(i), and their subtraction n-\{i) — 
n^{i) in a two-dimensional lattice Fermi gas, in which 
5 x 5 (= 25) square lattice is used, and U/t and V/t are 
taken —15 and 1, respectively. Although the system size 
is small, we find nj : ni ~ 2 : 1 inside the core, again. 
This result indicates that the penetration of un-paired 
excess atoms into the core region is a general feature in 
the presence of a lattice, irrespective of dimensionality of 
the system. 



VI. FINITE TEMPERATURE EFFECTS 

Finally, let us mention finite temperature effects on the 
approximate ratio 2:1 in the core phase. We find that 
there are two crucial temperatures in discussing the tem- 
perature effects on the present results. The first tempera- 
ture, which is much higher than the second temperature, 
is given by an energy scale in which the Hubbard model 
becomes sufficiently valid, i.e., the single-band Hubbard 
model is a good description. The experiments have al- 
ready succeeded in reaching a range below the first tem- 
perature. The second temperature is characterized by a 
lower energy range that the correlation effects have an 
important role, i.e., the present features become effec- 
tive. Thus, the problem is whether one can cool down 
below the second temperature or not. Let us estimate 
the second temperature. In the half-filling of 1-D lat- 
tice, 0.5t ~ E-p, where Ep is the Fermi energy and t is 
the the hopping matrix element in the Hubbard model, 
while the energy scale for the main correlation effect is 
t 2 /U as shown in Eq. ([3]) in this manuscript. Thus, if 
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one assumes U/t ~ 10, then the energy scale for the cor- 
relation ~ O.li. This means that the correlation effect 
becomes effective when the temperature is lower than 
0.2Ep. The temperature range T < 0.2i?F was achieved 
in the absence of the optical lattice. 

VII. SUMMARY 

We investigated ground state properties of two- 
component trapped Fermi gases with population imbal- 
ances loaded on optical lattices. Using the exact diago- 
nalization method and partly DMRG, we calculated the 
density profiles for 1-D and 2-D Hubbard model, as a 
function of the attractive interaction for various popula- 
tion imbalances. We showed that the population imbal- 
ance remains in the core phase, where the ratio of major 
to minor species, : ~ 2 : 1. This result is observed 
in the wide parameter range except for the weak attrac- 
tive interaction or the extremely strong trap one irrespec- 
tive of the dimensionality. This is different from the case 
in the absence of the lattice, where nj : nj. = 1 : 1 is 
observable in the core region. The character in the core 
phase on the optical lattice can be explained by the cor- 
relation effects on the lattice through the effective model. 

Note added — Recently, we became aware of Refs. [It], 
HH, which analyze the same model and give consistent 



results with ours. 
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